A COLLISION OF SUBCLUSTERS IN ABELL 754 
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ABSTRACT 

We present direct evidence of a collision of subclusters in the galaxy cluster Abell 754. Our comparison 
of new optical data and archival ROSAT PSPC X-ray data reveal three collision signatures predicted by 
n-body/hydrodynamical simulations of hierarchical cluster evolution. First, there is strong evidence of a 
non-hydrostatic process; neither of the two major clumps in the galaxy distribution lies on the off-center 
peak of the X-ray emission from the intracluster gas. Second, the peak of the X-ray emission is elongated 
perpendicular to the collision axis defined by the centroids of the two galaxy clumps. Third, there is evidence 
of comprcssion-hcatcd gas; one of A754's two X-ray temperature components (Henry & Briel 1995) is among 
the hottest observed in any cluster and hotter than that inferred from the velocity dispersion of the associated 
galaxy clump. These signatures are consistent with the qualitative features of simulations (Evrard 1990a, b) 
in which two subclusters have collided in the plane of the sky during the last ~Gyr. The detection of such 
collisions is crucial for understanding both the dynamics of individual clusters and the underlying cosmology. 
First, for systems Hke A754, estimating the cluster X-ray mass from assumptions of hydrostatic equihbrium 
and isothermality is incorrect and may produce the discrepancies sometimes found between X-ray masses 
and those derived from gravitational lens models (Babul & Miralda-Escude 1994). Second, the fraction of 
nearby clusters in which subclusters have collided in the last ~Gyr is especially sensitive to the mean mass 
density parameter Qq (cf. Richstone et al. 1992; Evrard et al. 1993; Lacey & Cole 1993). As we show for 
A754, we now have the means to identify recent collisions. With a large, well-defined cluster sample, it will 
be possible to place a new and powerful constraint on cosmological models. 
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c/j The complex dynamics of the gas, galaxies, and dark mat- 
C^ter in clusters of galaxies provides clues to cluster evolution 
and to the imdcrlying cosmology of the universe. In hi- 
erarchical models of the formation of large-scale structure, 
clusters evolve by accreting other clusters or smaller groups 
of galaxies from the field. Simulations show that when these 
systems collide, the effects can be observed for ~l/i^^ Gyr h 
is the Hubble constant in units of 100 km s~^ Mpc~^Evrard 
1990b). 

One key coUision signature is the existence of hot. X-ray 
emitting, intracluster gas that is not in hydrostatic equi- 
librium with the cluster potential. During a collision of 
two subclusters, the peaks of the X-ray emission and the 
galaxy surface density will not coincide because the gas, un- 
like the galaxies and dark matter, is collisional. Although 
the gross features of the optical and X-ray data agree in a 
few well-studied clusters (Fabricant et al. 1986, 1993), these 
samples lack the spatial resolution to test the critical scales 
< 0.5/i~^ Mpc where simulations predict differences. Today, 
multi-object spectroscopy and a new generation of X-ray 
telescopes enable us to compare the spatial distribution of 
hundreds of cluster galaxies with high-resolution X-ray sur- 
face brightness maps to test the assumption of hydrostatic 



equilibrium. 

Simulations show that the collision of the intracluster me- 
dia of the subclusters produces two additional signatures 
(Evrard 1990b, Roettiger et al. 1993). First, the gas as- 
sociated with a subcluster of galaxies is compressed in the 
direction of the collision, and if dense enough to be observed, 
appears as an elongated peak in the X-ray surface brightness 
distribution. Second, the collision significantly heats and 
may even shock some of the gas, creating an anomalously 
hot component in the spatially-resolved X-ray temperature 
map. The combination of high-sensitivity ROSAT PSPC 
data and multi-object spectroscopy allows us to test these 
predictions and to understand the history of recent collisions 
in clusters for the first time. 

We present the results from the first cluster we have ex- 
amined, Abell 754. There is evidence of a non-hydrostatic 
process; neither of the two major clumps in the galaxy distri- 
bution coincides with the off-center peak of the X-ray emis- 
sion. We also observe the other collision signatures. The 
X-ray peak is distorted in a direction perpendicular to the 
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axis on which the galaxy clumps lie, and the temperature 
of the X-ray gas surrounding one clump is a factor of two 
higher than that inferred from the clump's velocity disper- 
sion. 

We discuss the data briefly in §2 and the observed col- 
lision signatures in §3. We stress that deriving an X-ray 
cluster mass for a system that is neither isothermal nor in 
hydrostatic equilibrium is misleading. Most importantly, 
the detection of these signatures indicates that a coUision 
has occurred in A754 within the last Gyr. A deter- 

mination of the fraction of clusters with recent collisions, a 
more direct constraint of cosmological models than the frac- 
tion of clusters with substructure, is now within our grasp. 

2. THE DATA 

As part of a larger study of the internal dynamics of nearby 
clusters (Zabludoff 1995), we have obtained galaxy spectra 
and X-ray images for A754, a rich cluster at a mean rcdshift 
of 0.055 with an X-ray luminosity of Lx = 2.2 x lO'^^/j-^ 
erg s~^ (0.3-3.5 keV band; Fabricant et al. 1986). We used 
the Las Campanas multi-fiber spectrograph (Shectman et 
al. 1992) to acquire spectra of 481 galaxies with tub ^ 19 
that lie within 1° of the cluster center (Abell 1958). We 
determine radial velocities with a cross-correlation method 
(Tonry & Davis 1979) or by fitting emission lines. A com- 
parison with published data (Dressier & Shectman 1988) 
indicates that our velocities have la observational uncer- 
tainties of ~60 km s~^ and a zero point offset of < 60 km 
s^^. Cluster membership is assigned to 341 galaxies in the 
redshift range of 0.045 to 0.063 (limits corresponding to the 
cluster mean velocity 16400 ± (3 x the cluster velocity dis- 
persion of 900 km s~^)). The data and detailed analyses will 
be presented in later work (Zabludoff 1995); here we focus 
on comparing the positions of cluster members on the sky 
with the cluster X-ray surface brightness distribution. 

Wc acquired a 2322 sec exposure of A754 from the archives 
of the Position Sensitive Proportional Counter (PSPC) on 
the ROSAT satellite (Henry & Briel 1993). The X-ray emis- 
sion from A754 lies within the central 40 arcmin of the 2° 
PSPC field, where the point response function is approx- 
imately Gaussian, and the flatfield and background levels 
are approximately constant. Because A754 is detected in 
both the hard and soft PSPC energy bands, we use X-ray 
photons in the broad 0.1-2.4 keV band in our analysis to 
include all the signal. 

Plate 1 shows the smoothed X-ray map superimposed 
on the smoothed projected distribution of cluster galaxies, 
where each galaxy is equally weighted. We provide the de- 
tails, including the smoothing scales, in the Plate 1 caption. 
The plate is an unbiased map of the optical light in the 
cluster; weighting the galaxies by their luminosity does not 
significantly change the features. The peak of the X-ray 
emission is at (a,(5)2ooo ~ (09^' 09™ 22", -09° 40' 11"). "^c 
adopt an uncertainty in these coordinates of ~6 arcsec, the 
nominal attitude accuracy for telescope pointing given by 
the ROSAT mission description. By comparing the coor- 
dinates of the field centers and of the x-ray /radio galaxy 
(Harris et al. 1984) in the southwestern corner of the plate, 
we estimate that the errors in the alignment of the X-ray 
and optical contours are < 30 arcsec. 



Within the central Ih^^ Mpc, there is bimodal structure 
in the galaxy distribution. We assign a galaxy to either the 
southeast (SE) or northwest (NW) clump if the galaxy (1) 
lies within the clump's red optical contour in Plate 1 and (2) 
has a velocity within the limits corresponding to the clump 
mean velocity ± (3 x the clump velocity dispersion). Once 
memberships are assigned, we calculate the projected cen- 
troid of each clump. We estimate the uncertainty in the cen- 
troid position by throwing out each member, one at a time, 
and recalculating the centroid. The centroids vary during 
this procedure by < 10 arcsec. The SE galaxy clump lies 
at (a, (5)2000 ~ (09'' 09™ 26^ -09° 43' 49") and contains 41 
galaxies. The NW galaxy clump Hes at (a, 5)2000 ~ (09'' 08™ 
27^ -09° 37' 45") and contains 27 galaxies. The two clumps 
are separated by 15.6 arcmin (0.73/i~^ Mpc) on the sky. The 
mean velocities and velocity dispersions of the galaxy clumps 
aievNw = 16432 ±191 km s^^, aNW = 970^^34 km s^^ and 
vsE = 16320 ± 132 km s-\ asE = 830li?^ km s^^ (where 
V and (T are the bi-weight estimators of the line-of-sight lo- 
cation and scale; Beers et al. 1988). 

3. DISCUSSION 

A comparison of the X-ray and optical data in A754 yields 
three surprises. First, neither of the two major condensa- 
tions of galaxies corresponds to the most prominent feature 
of the highly asymmetric X-ray surface brightness map in 
Plate 1. The X-ray peak is offset from the SE galaxy clump 
by a projected separation of 3.8 arcmin (0.18/i^^ Mpc), 
which is much greater than the relative uncertainty in their 
positions (< 30 arcsec). Second, emission from the X-ray 
peak is elongated perpendicular to the axis joining the SE 
and NW optical clumps. Third, a comparison of the galaxy 
clumps and the distribution of X-ray gas temperatures re- 
veals that the two principal temperature components (Hen- 
riksen 1993; Henry & Briel 1995) are roughly aligned east- 
west along the axis connecting the optical clumps. Although 
the X-ray temperature of the eastern component (~5 keV), 
which surrounds the SE galaxy clump, is consistent with 
that derived from asE (5-2 keV: Tx = W'^ '^'^a'^'^-^ Edge & 
Stewart 1991), the temperature of the western component is 
~13 keV, much hotter than the 6.5 keV inferred from (Tnw 

The disturbed X-ray morphology is not due to the chance 
superposition of an unidentified X-ray source. The X-ray 
luminosity of the peak is ^20-25% that of the cluster. A 
single X-ray point source is ruled out by the elongation of 
the peak. The peak shape is not reproduced by a model 
of two point sources at any separation. Furthermore, no 
point sources are revealed in a deeper PSPC exposure of 
this field (Henry & Briel 1995). As for contamination from 
extended sources, there is no galaxy with < 19 within 
the brightest contour of the X-ray peak. The two closest 
galaxies, which are cluster members, lie near the southern 
and northern ends of the elongated 2nd- and 3rd-brightest 
contours, respectively. The southern galaxy, which has no 
emission lines and is classified SO/D (Dressier 1980), is visu- 
ally one of the brightest in the cluster {rriB = 16.3; Fabricant 
et al. 1986) and has a weak associated radio source (Harris 
et al. 1980). The northern galaxy (niB = 17.5) is typed Irr 
(Dressier 1980); its spectrum has [OH] ~ 40A > [OHI] ~ 
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lOA and so is consistent with a starburst. It is vcrv unlikelv 
that both galaxies have the X-ray luminosities (> lO^^h^"^ 
erg s~^) and X-ray radii (~100^~^ kpc) necessary to ac- 
count for the brightness and extent of the peak (Stocke et 
al. 1991; Fabbiano 1989). Visual inspection of the POSS 
plates finds no additional objects near the X-ray peak, ex- 
cept for an 8.9 magnitude F star (HD 78665) that makes 
a negligible contribution to the X-ray flux (Fabricant et al. 
1986). Thus, we conclude that the X-ray emission of the 
peak originates from the intracluster gas. 

What do the data for A754 imply about cluster evolu- 
tion? We claim that there has been a recent collision of 
two subclusters in A754 and now discuss three observations 
that have led to this conclusion. First, the displacement be- 
tween the galaxies, which we assume trace the dark matter, 
and the X-ray peak is consistent with that seen between the 
dark matter and the gas distributions during subcluster col- 
lisions in hierarchical clustering simulations (Evrard 1990b). 
During the collision, the models predict that the collisional 
X-ray gas is stripped from the free-streaming dark matter 
and galaxies. The gas is not in hydrostatic equihbrium with 
the dark matter potential for ~l/i~^ Gyr afterwards. The 
galaxies free-stream through the collision more easily than 
the gas, flying further than the gas from the point of impact 
after the clumps cross. It is suggestive that the SE galaxy 
clump and the X-ray peak lie in this conflguration and that 
their separation is consistent with the predicted displace- 
ments (< 0.5/i^^ Mpc). The NW galaxy clump does not 
have an obvious X-ray counterpart, but one is implied by the 
east-west eflipticity of the cluster's central X-ray isophotes 
and by the hot X-ray temperature component in the west. 

Second, the X-ray peak is elongated perpendicular to the 
collision direction. In simulations, this feature results from 
the compression of the gas after the intracluster media col- 
lide and persists for ~ lh~^ Gyr (Evrard 1990a, b; Roettiger 
et al. 1993). The high X-ray luminosity and relatively low 
X-ray temperature of the peak (compared to the hot western 
component) argue that it is a region of high gas density. 

Third, there are at least two temperature components in 
A754 (Henry & Bricl 1995), and the temperature of the west- 
ern component is anomalous. The ~12 keV gas in the west, 
which is much hotter than inferred from the velocity dis- 
persion of the galaxies in the NW clump, has been heated 
or even shocked beyond its equilibrium value, most likely 
by the compressive forces of the collision. Simulations pre- 
dict that this type of heating occurs after the clumps coUide 
(Roettiger et al. 1993), thereby supporting the picture in 
which the NW and SE galaxj^ clumps are not infalling for 
the first time, but have already crossed. 

The excellent agreement of the data with the qualitative 
appearance and scale of the models suggests that the galaxy 
clumps in A754 have crossed in the plane of the sky dur- 
ing the last ~l/i~^ Gyr. We can test this picture with 
a simple model of two subclusters on linear orbits (Beers, 
Geller, & Huchra 1982). For the input age of the uni- 
verse (8.28ft~^ Gyr; qq = 0.1), radial velocity difference be- 
tween the galaxy clumps (112 km s~^), clump projected 
separation {Rp = 0.73h^^ Mpc), and total system mass 
{Mtot = 4.8 X 10^'' Mq), there is a bound, outgoing solution 



of the motion equations in which the angle between the col- 
lision axis and the sky is 3.1 degrees, the clump separation is 
0.73h~^ Mpc (~ Rp), the relative clump velocity is 2071 km 
s~^, and the time elapsed since the clumps crossed is 0.2h~^ 
Gyr. The system mass Mtot, which is the virial mass esti- 
mated from the 57 cluster members within the circle whose 
origin is the centroid of the two clumps and whose diam- 
eter is Rp, is consistent with the sum of the virial masses 
of the two clumps. Changing Mfot by a factor of two does 
not affect the character of the solution. Because this model 
assumes that (1) the subclusters arc point masses, (2) the 
orbits are linear, (3) the influence of the perturbed intraclus- 
ter medium is negligible, and (4) the subcluster masses are 
constant, we use the model only to demonstrate the plausi- 
bility of our particular collision scenario. 

Plate 1 also provides insight into why cD galaxies in some 
clusters (Beers & Geller 1983; Zabludoff et al. 1990) do not 
lie on the projected, optical cluster center or at the bot- 
tom of the cluster potential well. The cD galaxy in A754 
(Dressier 1980) is within 1.2 arcmin of the centroid of the 
NW galaxy clump and has the same velocity as the clump 
mean velocity. The galaxy is therefore at rest in a local den- 
sity maximum, but is displaced relative to the projected, 
optical center of A754 as a whole. For subclusters on linear 
orbits in the plane of the sky, the cD's offset from the cluster 
center is maximized and its peculiar motion with respect to 
the cluster mean velocity minimized. We can imagine a less 
fortuitous configuration in which the collision axis is rotated 
90 degrees into the line-of-sight. In this case, both the bi- 
modality of the galaxy distribution and the displacement of 
the cD could only be detected in the distribution of galaxy 
radial velocities. The collision signatures would be nearly 
impossible to observe. Thus, any determination of the frac- 
tion of clusters with substructure or with recent subcluster 
collisions is a lower Hmit. 

A more complete understanding of the evolution of the 
subcluster collision requires n-body/liydrodynamical simu- 
lations that reproduce the specific optical and X-ray prop- 
erties of A754. New ASCA X-ray temperature maps with 
high spatial resolution may reveal the existence of a shock 
and the details of its structure. For now, we conclude that 
estimating cluster X-ray masses by assuming isothermality 
and hydrostatic equilibrium is perilous. In A754, the de- 
rived mass, assumed directly proportional to temperature, 
is overestimated by -^50% because the gas is heated beyond 
its equilibrium temperature (Henry & Briel 1993). This ef- 
fect may contribute to the discrepancies (Babul & Miralda- 
Escude 1994) sometimes found between cluster X-ray masses 
and the masses derived from gravitational lens models. 

How common are collisions and what are the implications 
for cosmology? In hierarchical models, the rate at which 
clusters accrete mass is especially sensitive to the mean mass 
density parameter flo (cf. Richstone et al. 1992; Evrard et 
al. 1993; Lacey & Cole 1993). Therefore, we can place 
stringent limits on Qq by determining the number of nearby 
clusters in which subclusters have recently collided. The ob- 
served fraction of clusters with substructure provides some 
estimate of the accretion frequency, but depends critically 
on (1) whether the apparent substructure is real or projected 
and (2) the time elapsed before the subclusters coUide and 
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virialize. This relaxation timcscalc in turn depends on the 
energy and eccentricity of the subcluster orbits and on the 
distribution of mass within the cluster. Because the relax- 
ation timescale is poorly constrained, it is not possible to 
know whether the substructure observed today reflects pre- 
dominantly recent or ancient accretion. This distinction is 
critical for the measurement of flo because the accretion rate 
is a strong function of lookback time. For example, the cur- 
rent fraction of clusters with substructure for an flo = 0.2 
model changes from 5% to 40% depending on whether the 
accretion occurred within the last ~2/i~^ or ^4:h~^ Gyr 
(Lacey & Cole 1993). As we show in A754, we can avoid 
most of these uncertainties by using the physics of both the 
gas and galaxies to directly identify clusters in which sub- 
clusters have collided in the last Gyr. 

We stress that A754 is unusual only because it has been 
examined on scales where the effects of subcluster collisions 
are most visible. Similar displacements between the galaxies 
and peak of the X-ray emission arc observed in a subcluster 
in Coma (Mellier et al. 1988; Briel et al. 1992) and in Virgo 
(Binggeh et al. 1993; Bohringer et al. 1994). It is notewor- 
thy that Coma and Virgo are the only two other clusters for 
which there are as many measured galaxy velocities as for 



A754. On-going collisions might also be the cause of incon- 
sistencies between the global kinematics of the gas and of 
the galaxies in Perseus (Gorenstein et al. 1978; Cowie, Hen- 
riksen, & Mushotzky 1987) and in Abell 2256 (Fabricant, 
Kent, & Kurtz 1989). A detailed comparison of the optical 
and X-ray properties of a large, well-defined sample of clus- 
ters will reveal the incidence of such colhsions and will place 
an improved constraint on cosmological models. 
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Plate Caption 

Plate 1 : The smoothed distribution of X-ray emitting gas (white contours) superimposed on the smoothed distribution of 
cluster galaxies in A754. The scale of the plot is 46.8/i~^ kpc arcmin"^ {qo = 0.1). Note that neither of the two major clumps 
of galaxies (bordered by red) hes on the extended peak of the X-ray emission. The distribution of 201 galaxies is smoothed 
with a Gaussian of 1.23 arcmin (~60/i^^ kpc) FWHM, a scale smaller than that of giant galaxy halos (Zaritsky & White 
1994). The faintest greyscale level corresponds to individual galaxies. The PSPC image is smoothed with a Gaussian of 30 
arcsec FWHM: the resulting map has a resolution of -^40 arcsec (~30/i~^ kpc) determined from the FWHM of the X-ray 
point source ~8 arcmin north of the X-ray peak. The background level is 7.0 ±0.1 x 10"" counts s~^ arcmin"^. The lowest 
contour level is four times the background. The highest level plotted is 0.04 counts s^^ arcmin^^ (57 x the background). The 
difference between each contour is 2.9 x 10^^ counts s^^ arcmin^^. The bright X-ray source in the southwest is associated 
with a radio galaxy in a group of galaxies at the same redshift as A754 (Harris, Constain & Dewdney 1984). 



